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ABSTRACT 

The early gas-dust solar nebula is considered: the gasdynamic theory is used 
to study the gravitational Jeans-type instability in its protoplanetary disk. The 
implications for the origin of the solar system are discussed. It is shown that a 
collective process, forming the basis of the gravitational instability hypothesis, 
solves with surprising simplicity the two main problems of the dynamical charac- 
teristics of the system, which are associated with its observed spacing and orbital 
momentum distribution. 

Subject headings: planetary systems: formation-solar system: formation 
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1. Introduction 

Many young stars are surrounded by gas-dust disks (Bodenheimer & Lin 2002). 
Planetary formation is thought to start with inelastically colliding gaseous and dust particles 
settling to the central plane of a disk to form a thin and relatively dense layer around 
the plane. During the early evolution of this disk it is believed that the dust particles 
coagulate into kilometer-sized rocky asteroids- "planetesimals" (~ 10 10 such bodies) owing 
to the gravitational instability (Goldreich & Ward 1973) and/or to the collisional sticking 
(Beckwith et al. 1990). Of these processes, dust particle settling can now be observable. 

We suggest that all planets of the solar system were created by disk instability. That 
is, as a result of local gravitational instability, on attaining a certain critical thickness 
(and density, respectively), small in comparison with the outer radius of the system R, 
the circumsolar gas-dust disk disintegrated into a large number of separate protoplanets. 
Following Boss (2003), this hypothesis envisions coagulation and settling of dust grains 
within the protoplanets to form rock and ice cores. A protoplanet accreted a gas 
subsequently from the solar nebula after accumulating a solid core of ~ lM e , followed 
by the loss of the light elements of the terrestrial planets through the thermal emission of 
the sun. The advantages of the disk instability model are that (1) the instability process 
itself is quite fast, and could form planets in 10 3 — 10 4 yr (Boss 2002) and (2) in unstable, 
nonaxisymmetric disks differential rotation can simultaneously transfer angular momentum 
outward and mass inward through gravitational torques. The work described here has 
precedents in earlier studies of gravity disturbances in galactic disks and the Saturnian 
ring disk (e.g., Shu 1970; Lynden-Bell & Kalnajs 1972; Griv, Yuan & Gedalin 1999; Griv, 
Gedalin & Yuan 2003). 
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2. Dispersion relation 

Let us consider the dynamics of the gaseous component in the presence of the collective 
self-gravitational field. A Langrangian description of the motion of a fluid element under 
the influence of a spiral field is used, looking for time-dependent waves which propagate in 
a differentially rotating, two-dimensional disk. The approximation of an infinitesimally thin 
disk is a valid approximation if one considers perturbations with a radial wavelength that 
is greater h, the typical disk thickness (Toomre 1964; Shu 1970; Genkin & Safronov 1975; 
Safronov 1980). 

The time dependent surface density H(r,t) is splited up into a basic and a developing 
(perturbation) part, £ = £ (r) + £i(r, t) and |Ei/E | <C 1, where r, ip, z are the cylindrical 
coordinates and the axis of the disk rotation is taken oriented along the z-axis. The 
gravitational potential of the disk K(r, t) is also of this form. These quantities £ and K are 
then substituted into the equations of motion of a fluid element, the continuity equation, 
the Poisson equation, and the second order terms of the order of Sf, may be neglected 
with respect to the first order terms. The resultant equations of motion are cyclic in the 
variables t and </?, and hence by applying the local WKB method one may seek solutions in 
the form of normal modes by expanding any perturbation 

Ei(f,*),Ni(f,f) = 5Z,5Ke lkrr+imip - wt + c.c, (1) 

where 5T, and 5ti are the real amplitudes, which are constant in space and time, k r (r) 
is the real radial wavenumber, m is the nonnegative (integer) azimuthal mode number, 
uj = Sfta; + i^scu is the complex frequency of excited waves, and c. c. means the complex 
conjugate. The solution in such a form represents a spiral plane wave with m arms. The 
imaginary part of uj corresponds to a growth (Quj > 0) or decay (Quj < 0) of the components 
in time, Si and oc exp(Qwt), and the real part to a rotation with constant angular 
velocity Q p = 3t,oj/m. Thus, when > 0, the medium transfers its energy to the growing 
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wave and oscillation buildup occurs. 

It is important to note that in the WKB method, the radial wavenumber is presumed 
to be of the form 

kr(r)=AV(r), (2) 

where A is a large parameter and \l/(r) is a smooth, slowly varying function of the radial 
distance r, i.e., dln/c r /dlnr = 0(1), and \k r \r ^> 1. 

Paralleling the analysis leading to equation (34) in Griv et al. (1999), it is 
straightforward to show that 

3E 30 2 + w> 2 2fim<91nE \ 

— K V r oo* r or J 

where Ei(£ — > — oo) = 0, so by considering only growing perturbations we neglected the 
effects of the initial conditions, uj* = uj — mfl is the Doppler-shifted (in a rotating reference 
frame) wavefrequency, f2(r) is the angular velocity of differential rotation at the distance r 
from the center, and k ~ is the epicyclic frequency. In equation above, 3 = Ki + Pi/Eo, 
P\ is the perturbed gaseous pressure, and c = (dP/ '(9E) 1 / 2 is the sound velocity. In 
equation (3) only the most important low-frequency (|u; 2 | < k 2 ) perturbations developing 
in the plane z = between the inner and outer Lindblad resonances are considered (Griv et 
al. 1999, 2003). Equating the perturbed density Ei [eq. (3)] to the perturbed density given 
by the asymptotic (k 2 3> m 2 /r 2 ) solution of the Poisson equation (Griv et al. 1999), the 
Lin-Shu-type dispersion relation is obtained 

O m 

W*l,2 ~ ±p|wj| - 27tGE — —rj- , (4) 

OJj T | K | Li 

where p — 1 for gravity-stable perturbations with u 2 ~ u 2 > 0, p = i for gravity-unstable 
perturbations with < 0, L = (9 In E /<9r) _1 is the radial scale of spatial inhomogeneity, 
\kL\ 3> 1, and the term involving L~ l is the small correction. Also, 

tu 2 = K 2 -2nGi; (k 2 J\k\) + k 2 c 2 (5) 
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is the squared Jeans frequency, k = \/k 2 + m 2 /r 2 is the total wavenumber, k 2 = 

k 2 {l + \{2Vt/ k) 2 — 1] sin 2 -?/'} is the squared effective wavenumber, and ip = arctan(m/r/c r ) 

is the perturbation pitch angle. 

Equation (4) determines the spectrum of oscillations. In the gravity-unstable case, the 
equilibrium parameters of the disk and the azimuthal modeazimuthal mode number m (= 
number of spiral arms) determine the spiral pattern speed of Jeans-unstable perturbations 
(in a rotating frame): 

Q p = $tuj*/m pa 27rG£ -^- ryry , (6) 

I ^ j I ^ I k I L' 

where 27rG , S |A;| ~ fi 2 , |cj 2 | ~ fi 2 , rk 2 \L\ > 1, and, therefore, fi p ~ Vt/rk 2 L < Q. Thus, the 
typical pattern speeds of spiral structures in Jeans-unstable, oj 2 < 0, disks are only a small 
fraction of some average angular velocity fi av . Because fl p does not depend on m, each 
Fourier component of a perturbation in an inhomogeneous system will rotate with the same 
constant angular velocity. The theory states that in homogeneous (|L| — > oo) disks fi p = 0. 

The disk is Jeans-unstable to both axisymmetric (radial) and nonaxisymmetric (spiral) 
perturbations if c < c T , where c T = ttGT^q/k is the Safronov-Toomre (Safronov 1960, 1980; 
Toomre 1964) critical velocity dispersion to suppress the instability of axisymmetric (ip = 0) 
perturbations. Thus, if the disk is thin, c rfl, and dynamically cold, c < ct, then such 
a model will be gravitationally unstable, and it should almost instanteneously (see below 
for a time estimate) taken on the form of a cartwheel. The instability, which is algebraic 
in nature, is driven by a strong nonresonant interaction of the gravity fluctuations (e.g., 
those produced by a spontaneous perturbation and/or a satellite system) with the bulk of 
the particle population, and the dynamics of Jeans perturbations can be characterized as a 
nonresonant interaction, that is, in equation (3), uj^ — Ik^ 0, where / = 0, ±1. 

A very important feature of the instability under consideration is the fact that it is 
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almost aperiodic (|$Rw*/9iu;*| Cl). The growth rate of the instability is relatively high: 

Stu* « ^2%GZo(k 2 J\k\) (7) 

and in general 3w* ~ Q, that is, the instability develops rapidly on a dynamical time 
scale (on a time of 3 — 4 disk rotations, or about 10 4 yr in the early solar nebula). 
From equation (5), the growth rate of the instability has a maximum at the wavelenght 
A cr it ~ 2c 2 /GS . At the boundary of instability, that is, c ~ or, A cr i t ~ 27r 2 GS //t 2 ~ 2nh. 
It means that of all harmonics of initial gravity perturbation, one perturbation with 
A crit ~ 2nh, with the associated number of spiral arms m, and with the pitch angle ip 
will be formed asymptotically in time af a single rotation 5 x 10 9 yr ago). For the 
parameters of the solar nebula (R ~ 300 AU, k = 2n/T orh ~ 10 -10 s -1 , and the total 
mass of the disk M d ~ 0.1 M Q ), one obtains the typical mass of the core of a protoplanet 
M c ~ 10 "U. ~ M e . 



3. Spacing of the planets 

There exists the empirical Titius-Bode (TB) rule which gives the mean orbital distances 
of the planets and which can be written in the Blagg-Richardson formulation as 

r n = r A\ (8) 

where r n is the distance of the nth planet from the Sun (in AU), n = 1 for Mercury, 2 
for Venus, . . ., and 9 for Neptune, A = 1.73 is the mean ratio between two consecutive 
planetary distances, and r ~ 0.21. Also, one cannot overlook the fact that many of the 
regularities which are found in the planetary system are also to be seen in the regular 
satellite systems of Jupiter, Saturn, and Uranus, e.g., the spacing of the regular satellites 
is a variation of the TB rule (Fig. 1). This suggests that the same cosmogonic process 
must have been responsible for the origin of both types of systems. Lynch (2003) has 
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already argued that it is not possible to conclude unequivocally that laws of TB type are, 
or are not, significant. Therefore, the possibility of a physical explanation for the observed 
distributions remains open. 1 
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Fig. 1. — Relation between distances of planets (satellites) from the Sun (planets) r and 
their numbers n. Data observed are represented by circles: (a) the solar system, (b) the 
satellite system of Jupiter, r = 249.679 and A = 1.649, (c) the satellite system of Saturn, 
r = 92.416 and A = 1.503, and (d) the satellite system of Uranus, r = 89.737 and A = 1.46. 
The crosses represent the TB rule, equation (8). 

Equation (8) can be rewritten: 

(2tt/ In 1.73) ln(r n /0.21) = 2™ . (9) 

Next, the surface density of the disk may be represented in the form of the sum of the 
equilibrium surface density S (r) and the perturbed surface density 



£i(r) =SE(r)e Jut cos [11.5 \n(r n /0.21) + rrvp] , 



(10) 



interestingly, the mean orbital distance to the recently discovered classical Edgeworth- 
Kuiper belt objects, r 46 AU, is in fair agreement with that given by the TB rule for the 
solar system's 9th planet, r 10 « 50 AU. 
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where <5£(r) is the amplitude varying slowly with radius, and [11.5 ln(r n /0.21) + mip] 
represents the phase varying rapidly with radius, 

\k r \r = 11.5|(d/dr) ln(r n /0.21)| r > 1 . 



Equation (9) and the condition 5E(r) > on the initial phase imply that the maximum 
values of the perturbed density in equation (10) coincide with the positions of all the 
planets (Fig. 2a). 




5 



r, radius (AU) 

Fig. 2. — (a) Dependence of the perturbed surface density of the protoplanetary disk £i(r) 
(arbitrary units) on the radius r, equation (10). The maxima of the perturbed density 
coincide with locations of all the planets, (b) Spiral density waves with m = 1 arm [eq. (10)] 
in the (r, </?)-plane, (c) density waves with m = 2 arms, and (d) density waves with m = 3 
arms. The filled circles represent the maxima of the perturbed density (protoplanets) of 
Jeans waves, which are unstable to both axisymmetric and nonaxisymmetric perturbations. 

Interestingly, and this is the central part of our theory, the TB rule [eq. (10)] satisfies 
the conditions of the WKB wave with the effective TB radial wavenumber 

11.5 , . 

keft = , (11) 

dln£; cfr /dlnr = 0(1), and k cS r > 1 [cf. eq. (2)]. 2 

2 Polyachenko (Polyachenko & Fridman 1972) has already been considered this analogy 
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Thus, if the space dependence of the perturbed surface density of the protoplanetary 
disk in the (r, </?)-plane has the form of equation (10) with > 0, the maxima of both 
radially and azimuthally unstable gravity perturbations are located in places of the solar 
system's planets (Figs 2b, c, d). Let us define conditions under which the density maxima 
are localized on planetary orbits. If the disk is inhomogeneous with respect to equilibrium 
parameters, the wavelength of a perturbation with a maximum growth rate A cr jt will be a 
function of the radius r. From the above, the wavelength A crit ~ 47r 2 G , £ //t 2 , corresponding 
to the minimum on the dispersion curve (4) (see also Griv et al. 2003, Fig. 1 therein). On 
the the other hand, the wavelength is A c g = 2ir/k c fi. Comparing A cr i t with A e ff , we see that 
in the case where the disk density is dependent on radius according to the law 

£ (r) « 0.0138G" Vr , (12) 

the maxima of time-increasing, both radially and azimuthally Jeans-unstable density 
perturbations are arranged in it according the TB rule given by equation (8). 

The last condition may be fulfilled in Keplerian disks, k oc r _3//2 , only if S oc r~ 2 . 
Interestingly, Tomley et al. (1991) have used almost the same law S oc r~ 7 / 4 as initial 
profile for simulation of a disk surrounding the central star. The reason for using such 
a law comes from a particular model of protostellar cloud collapse Tomley et al. used. 
It was obtained that this initial model did not get much subsequent evolution in the 
simulations although it was nicely gravitationally unstable. Based on hydrostatic models, 
the radial density distributions in circumstellar disks around Herbig Ae/Be and T Tauri 
stars have been proposed to be in the range S oc r _(1 - 9 " 2 - 4 ). Detailed modeling of the 

in his investigation of the possibility of the explanation of the law of planetary distances by 
the gravitational instability in sufficiently flat systems, but evidently without success. In 
particular, Polyachenko studied only axisymmetric m = perturbations, which do not carry 
angular momentum (see an explanation below in §4). 
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NIR-to-millimeter appearance of several spatially resolved T Tauri disks has confirmed these 
predictions. It has been stated that optically thick young disks around those stars with 
spatial structures are dominated by gravitation and gasdynamics. See, e.g., Stapelfeldt et al. 
(1998), Chiang & Goldreich (1999), and Wolf, Padgett & Stapelfeldt (2003) for a discussion. 
Fits of models to observed spectral energy distributions of protostellar disks typically give 
S oc r~ 3 / 2 (Bodenheimer & Lin 2002). Also, a standard reference model of a disk, known 
as the "minimum mass solar nebula," reconstructed from the distribution of mass in the 
planets of the solar system and assuming solar composition and no migration of planets, 
gives S oc r~ 3 / 2 . The latter is close to the r~ 2 distribution advocated above. Clearly, 
given the observational and analytical uncertainties, the two distributions, S oc r~ 3 / 2 and 
S oc r~ 2 , are not necessarily inconsistemt with each other. For instance, the inclusion of 
the disk's self-gravity in addition to the gravitational field of the sun will reduce the value of 
the exponent n in the required density-radius relation S oc r~ n . In turn, both optical and 
near-infrared observations of pre-main-sequence stars of intermediate mass have revealed 
the spiral structure, and thus presumably the Jeans instability, in the circumstellar disks 
(Grady et al. 1999; Clampin et al. 2003; Fukagawa et al. 2004). 

One concludes, therefore, that if the surface density of a protoplanetary disk falls 
according to the law given by equation (12), the increasing maxima of density perturbations 
of a Safronov-Toomre-unstable disk (c < ct) are located between the Lindblad resonances 
in places of the planets (Fig. 2). We believe to have obtained a theoretical interpretation 
of the TB rule: the distance between planets is the wavelength of the most Jeans-unstable 
perturbations at the given point of the protoplanetary disk. 

By using equation (12), it is easy to find that the disk mass between 0.3 AU and 30 AU 
is about 0.4 M . This means that in the present planets there is contained not more than 
about 0.5% of the mass of the protoplanet cloud. Almost certainly, a part of the initial 
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mass of the planets was blown away due to intensive corpuscular emission of the early sun. 



4. Orbital momentum distribution 

We next turn to the question of how to account for the concentration of angular 
momentum in the planets and of mass in the sun. The torque exerted by the gravity 
perturbations on the disk is T = — J J d 2 r(r x VNi)£i or 

r—r^r^r,^i^. (13) 



j ri Jo d(p' 
The points r\ and r 2 in which a>* ± k = are called the points of inner and outer Lindblad 
resonances. They play an important role in the theory: the solution of spiral type (1) 
rapidly oscillating in the radial direction lies between r 1 and r 2 . Outside the resonances, 
r < r\ and r > r 2 , the solution decreases exponentially. A special analysis of the solution 
near corotation (u>* = 0) and Lindblad resonances is required. Resonances of a higher order, 
cj* ± Ik = and |/| = 2, 3, • • •, are dynamically of less importance (Shu 1970). To emphasize 
it again, the present analysis is restricted to consideration of only the principal part of 
a disk between the Lindblad resonances. Investigation of the wave-particle interaction at 
spatially limited resonances has been done by Lynden-Bell & Kalnajs (1972), Goldreich & 
Tremaine (1978, 1980), and Griv, Gedalin, Eichler & Yuan (2000). 

Using equation (3), from equation (13) one finds 

where Qw* > 0, is the complex conjugate potential, and the values of Ki, K^, S , £1 are 
evaluated at r = r\. Three physical conclusions can be deduced from equation (14). First, 
the distribution of the angular momentum of a disk will change under the action of only 
the nonaxisymmetric forces oc m. The latter is obvious: axially symmetrical motions of a 
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system, studied by Polyachenko, produce no gravitational couples between the inner parts 
and the outer parts. Second, the distribution of the angular momentum will change upon 
time only under the action of growing, i.e., Jeans-unstable perturbations (Qw* > 0). 3 Third, 
T < 0: the spiral perturbations remove angular momentum from the disk. This takes place 
in the main part of the disk between the Lindblad resonances where spiral density waves 
are self-excited via a nonresonant wave- "fluid" interaction. Further there is absorption 
of angular momentum by particles that resonate with the wave (Lynden-Bell & Kalnajs 
1972). As a result, the bulk of angular momentum is transferred outward (and a mass 
transported inward, correspondingly). In turn a small group of resonate particles moves 
outward taking almost all angular momentum. 4 These processes lead to the core-dominated 
mass density profile in the protoplanetary disk, together with the buildup of an extended, 
rapidly rotating outer envelope. We speculate that a large portion of the initial mass of the 
nebula was transported toward the sun. 

Let us evaluate the gravitational torque for a realistic model of the protoplanetary 
disk. In accordance with the theory developed above, the fastest growing spiral mode with 
m > 1, k* — k crit , and So;* ~ Q is considered. Taking into account that 87rm 2 K!N^ ~ Kg 
(an astrophysicist might well consider a perturbation with Ki/K of 1/10 or even 1/3 to be 
quite small) and K ~ r 2 Q 2 , where K is the basic potential, from equation (14) one obtains 
|T] ~ X r 4 f2 2 . The angular momentum of the disk C ~ £ r 4 f2. Then the characteristic 

3 In the opposite limiting case of slow growth (So;* — > 0), absorption and emission of 
angular momentum are confined only to resonate particles (e.g., Lynden-Bell & Kalnajs 
1972). The treatment of resonances is beyond the scope of the present analysis. 

4 Lynden-Bell & Kalnajs (1972) have proved that in good conformity with TV-body simu- 
lations the gravitational torques can only communicate angular momentum outward if the 
spirals trails. 
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time of the angular momentum redistribution is t ~ C/T ~ Q^ 1 . Thus, already in the 
first 3 — 4 disk revolutions, in say about 10 4 yr, the gas-dust protoplanetary disk sees its 
almost all angular momentum transferred outward and mass inward. We conclude that the 
Jeans instability studied here can give rise to torques that can help to clear the nebula on a 
time scale of > 1 Myr, in accord with astronomical requirements. In addition, the analysis 
is found to imply the existence of a new planet (or another Kuiper-type belt) at a mean 
distance from the sun of r n = 0.21 x 1.73 11 w 87 AU. 
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